A number conserving theory for topologically protected degeneracy in 

one-dimensional fermions 
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Semiconducting nanowires in proximity to superconductors are among promising candidates to 
search for Majorana fermions and topologically protected degeneracies which may ultimately be used 
as building blocks for topological quantum computers. The prediction of neutral Majorana fermions 
in the proximity-induced superconducting systems ignores number-conservation and thus leaves open 
the conceptual question of how a topological degeneracy that is robust to all local perturbations 
arises in a number-conserving system. In this work, we study how local attractive interactions 
generate a topological ground-state near-degeneracy in a quasi one-dimensional superfluid using 
bosonization of the fermions. The local attractive interactions opens a topological quasiparticle gap 
in the odd channel wires (with more than one channel) with end Majorana modes associated with 
a topological near-degeneracy. We explicitly study the robustness of the topological degeneracy 
to local perturbations and find that such local perturbations result in quantum phase slips which 
split of the topological degeneracy by an amount that does not decrease exponentially with the 
length of the wire, but still decreases rapidly if the number of channels is large. Therefore a bulk 
superconductor with a large number of channels is crucial for true topological degeneracy. 

PACS numbers: 



I. INTRODUCTION 

In the past few years, topological superconductors 
have become promising candidates for realizing Majorana 
fermions (MFs)i"— , which are of interest both for funda- 
mental reasons as a new type of particle with non-Abelian 
statistics and also their potential application in topo- 
logical quantum computation (TQC)^^"— . One of the 
simplest such topological superconducting (TS) systems 
supporting MFs consists of a semiconductor nanowire in 
a magnetic field together with spin-orbit coupling and 
superconductivitj^ii^ii^ through the appropriate tuning 
of the semiconductor chemical potential or equivalently, 
the carrier density. The superconducting pair-potential 
is induced in the nanowire by the proximity effect from 
an s-wave superconductor in contact with the nanowire. 
It has been shown that such a nanowire can be driven 
into a TS phase with MF end modes^^ii^^i^. 

The MFs in the TS phase of the semiconducting 
nanowire are localized at the ends of the wire and can be 
detected as zero-bias conductance peaks in the tunneling 
spectrum at the end o^°'^^i^° . The end MFs are associated 
with self-adjoint operators 71 and 72 (7I = 7i) and in this 
sense MFs are their own anti-particle. MFs are of inter- 
est for TQC because despite having no internal degrees of 
freedom individually, a pair of MFs has two distinct pos- 
sible states (fusion channels). These states are energeti- 
cally degenerate to a degree exponential in the separation 
of the MFs, and correspond to states where the combined 
fermionic state = ^^^^ is occupied (particle fusion 
channel) or empty (vacuum fusion channel) . These chan- 
nels correspond to even and odd total number of fermions 




Zeeman splitting 



FIG. 1: (Color online) (a) Combination of Zeeman splitting, 
spin-orbit coupling and weak attractive s-wave Feshbach reso- 
nance can induce end MFs in multi-channel nanowire (shown 
as green rectangle) . Spin-orbit coupling and Zeeman splitting 
ensures almost helical fermions at the fermi surface (shown as 
discs with arrows representing spin). Attractive s-wave inter- 
action leads to the scattering of pairs of fermions with zero 
momentum (Cooper pairs) from one channel (represented by 
blue color) to another (represented in red), (b) Schematic 
band-structure of multi-channel nanowire with all degenera- 
cies lifted by Zeeman splitting. Filling corresponding to odd 
number of bands is in the topological phase. 



in the system. The MFs in a system of nanowires have a 
non-Abelian statistics under exchange of positions of the 
MFs^ir— . The result of exchange of MFs in the nanowire 
system, which can be determined by measuring the chan- 
nel (particle or hole) each pair of MFs fuses into, de- 
pends on the order in which the exchanges are performed. 
This non-Abelian statistics forms the basis of topologi- 
cally protected operations using MFs. Most studies of 
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MFs in TS have so far been restricted to non number- 
conserving cases either via mean-field BCS studies or 
for systems where the superconductivity is induced by 
the proximity effect. Within these formahsms the MF is 
necessarily neutral. However tunneling of electrons into 
MF modes results in a transition of the nanowire between 
the quasiparticle and vacuum fusion channels which have 
different total electron numbers. Therefore a proper de- 
scription of tunneling through a neutral MF mode must 
involve charge fluctuations where the initially localized 
electronic charge delocalizes over the whole wire. A re- 
cent study by Leggett^^ discusses the issue of manifestly 
including number conservation and constructs a number- 
conserving version of a p-wave topological superconduct- 
ing wire HamiltoniarJ^ with topological degeneracy for a 
particular set of parameters. 

In this paper, we analyze the consequences of num- 
ber conservation by comparing several models in which 
we have one or more coupled nanowires in the TS state 
where superconductivity is generated by intrinsic at- 
tractive interactions. The calculations in this paper 
are performed within the bosonization technique2Ii2§,, so 
that any recourse to mean-field BCS-type approxima- 
tions that break number conservation is avoided. Each 
wire is supposed to contain an odd number Nch of non- 
degenerate, spin-orbit split, transverse channels. In our 
calculations, we assume that the wires may have an in- 
trinsic superconductivity, induced by attractive two-body 
interactions of finite range. We will also comment on 
the comparison of this scenario to the case where su- 
perconductivity is induced by the proximity effect from 
a neighboring bulk superconductor. As such our work 
will also apply to the one-dimensional topological insu- 
lator nanowire system which also has an odd number of 
conducting channels.— Our focus will be on enumerat- 
ing the low energy states of the wire network, associ- 
ated with MF modes at the ends or at junctions between 
wires, and on estimating the residual energy splittings be- 
tween the low energy states which occur when the wires 
are long but not infinite in length. In particular, we 
will want to distinguish between energy splittings that 
fall off exponentially in the length of the wire segments, 
and those which have a weaker dependence on length. 
Our analysis will use bosonization techniques appropri- 
ate for one-dimensional fermions, together with renor- 
malization group and scaling analyses where needed.— 
The bosonization technique has already been applied in 
the context of TS nanowires to understand the effect 
of interactions^ in non-number conserving proximity- 
induced superconducting systemSi ^^'^" 

We note that a few-channel nanowire^ with lo- 
cal attractive interactions can be potentially be real- 
ized in cold-atom systems with quasi-one-dimensional 
confinement, together with Zeeman spitting, spin-orbit 
coupling22, and attractive s-wave Feshbach interactions^. 
However, the primary motivation for this paper is to ex- 
plore the conceptual issues related to conservation of par- 
ticle number^^ and how topological degeneracy may be 



described in such a framework. 

In the case of proximity-induced superconductivity, it 
has typically been assumed that the bulk superconduc- 
tor is infinitely large, so that the phase of the supercon- 
ducting order parameter may be perfectly well defined, 
and one does not need to worry about conservation of 
total charge. Moreover, it has generally been assumed 
that bulk superconductors attached to different wires are 
strongly coupled to each other, so that they all have the 
same superconducting phase, and one does not have to 
worry about effects of phase slips. In our discussions be- 
low, we shall relax this last assumption and comment 
on the situation when phase slips do occur in the con- 
nections between different wires, either due to thermal 
activation or due to quantum-mechanical tunneling. For 
the case of nanowires with intrinsic superconductivity, we 
shall primarily be interested in networks of several wires, 
which are coupled to each other through tunnel barriers 
over large portions of their lengths. Results will depend 
on the number of wires involved and the ways in which 
they are coupled. 



II. OUTLINE AND RESULTS 

We first consider the case of single isolated wire. To 
clarify the role of number conservation, we start in 
Sec.|TTT]by discussing the differences between the degen- 
eracies of proximity-induced TS wires and intrinsic TS 
wires. Following this discussion, it becomes clear that 
for an isolated wire, the size of the energy gap will not 
decrease exponentially with the length L of the wire, but 
will generally be inversely proportional to L. This is a di- 
rect consequence of particle conservation in the isolated 
wire. Similar splitting has been predicted for MFs on 
superconducting grains.=34 

In order to find a situation where the tunneling gap 
and the lowest particle-hole excitation fall off faster than 
we must consider at least two coupled wires. Before 
discussing multi-wire systems, we start in Sec. llVl bv re- 
viewing how to determine the properties of weakly inter- 
acting multi-channel wires using bosonization. In SecfVl 
by computing the scaling dimensions of the Cooper-pair 
scattering terms, we find that for an infinite wire, with 
Nch > 1, for weak interactions which are attractive in all 
channels (see more precise definition below) the system 
develops an energy gap for single electron tunneling at 
points away from the ends of the wire, along with quasi- 
long-range power-law correlations in the superconducting 
order parameter As we will discuss in Sec. |Vl despite 
the absence of a global superconducting phase, bosoniza- 
tion allows us to assign a local superconducting phase 
29a{x) to each of the Nch channels. In this phase, the 
values of 6a (x) are locked to each other up to a multiple 
of TT. Thus, the phases of each channel relative to the 
local common superconducting phase forms a Z2 order 
similar to the Ising model. In fact, at this level, bosoniza- 
tion maps the multi-channel wire into Nch copies of spin- 



3 



less p-wave TS wireai* coupled to a global phase field, 
each of which can be mapped to a transverse field Ising 
model by a non-local Jordan- Wigner transformation.— 
However, this Z2 phase representation does not mani- 
festly respect the particle number conservation symme- 
try of the intrinsically number conserving interacting TS 
Hamiltonian. This is remedied in Sec. IVII bv identifying 
the appropriate total number subspace which is preserved 
by the Hamiltonian. The restriction of the Hilbert space 
reduces the topological degeneracy from 2^'='' to 2^'=''"^. 

The degeneracy associated with Zi order parameters 
for each channel is topologically protected only if it is 
robust to all local perturbations in the fermionic repre- 
sentation. Bosonization is a non-local transformation, 
similar to the Jordan- Wigner transformation, and there- 
fore, the local Zi order associated with each channel in 
the bosonic representation is a non-local and potentially 
topological order parameter in the original fermionic rep- 
resentation. However, the topological robustness to local 
fermionic perturbations is something that must be explic- 
itly checked. Since local fermionic tunneling operators 
are gapped in the middle of the wire, the ends of the wires 
are expected to be more susceptible to fermionic pertur- 
bations. By considering all such local fermion terms at 
the ends of the wires, in the bosonic representation in 
Sec. IVHl we show explicitly in Sec. IVHIl that the topo- 
logical degeneracy is robust to perturbations at the ends 
of the wire only for wires with an odd number of chan- 
nels Nch- In the absence of the number conservation 
constraint, the end scattering reduces the degeneracy of 
an isolated wire from 2^"^ to 2^r_h[™-odT) ^ This is similar 
to the familiar splitting of pairs of MFs in the proximity- 
induced case.^^ However, since in this paper we are inter- 
ested in finding exponential in system size splittings of 
topological degeneracies, we derive this result entirely in 
the bosonic representation without recourse to refermion- 
ization. Thus, topological degeneracy can only survive 
for Nch odd. Of course, for the number conserving case, 
the degeneracy is completely eliminated and one must 
use a bundle of wires to obtain a degeneracy. 

The robustness to scattering in the middle of the is 
a more subtle issue that is clarified later. For a wire 
with ends, the topological degeneracy manifests itself as 
states with energies inside the bulk gap, which can be ac- 
cessed by tunneling of an electron in the vicinity of either 
end. The amplitude for tunneling to these states falls off 
exponentially as the tunneling point moves away from 
either end. For a wire with an odd number of channels, 
the energy gap for tunneling near the ends of a wire will 
decrease as the wire length increases, and the relevant 
states will be reminiscent of the Majorana end states for 
a wire with proximity-induced superconductivity. 

As will become clear from both the general consid- 
erations in Sec. IHII and also the detailed bosonization 
treatment in the following sections, in order to find a sit- 
uation where the tunneling gap and the lowest particle- 
hole excitation fall off faster that we must consider 
at least two coupled wires as discussed in Sec. IIXI The 



simplest case is to consider parallel wires of length 
L, each with an odd number of channels, which are cou- 
pled to each other only over the middle third of their 
length. (Fig. [^b)) If the tunnel coupling between the 
wires is sufficiently strong, we expect that the number of 
electrons in each separate wire will fluctuate by a large 
amount, (with root-mean-square fluctuation proportional 
to L~^l'^). Then the fermion parity in any particular wire 
will have little effect on the energy, as long as the total 
fermion parity is fixed. Thus the energy difference be- 
tween the lowest state in which say wires 1 and 2 have 
both have even parity, or they both have odd parity, 
could be quite small. Indeed, we find that under ideal 
conditions, where the tunneling between the wires turns 
on adiabatically as a function of the distance along the 
wires, the energy splitting between lowest energy states 
of the system can be exponentially small in L. This is 
shown in Sec. |X] by analyzing the role of instanton tun- 
neling terms where the global phase of an entire channel 
can jump by tt from one energy minimum to another. 
The requirement of strong coupling between two wires 
is achieved when the Josephson coupling energy between 
the wires is large compared to the charging energy of ei- 
ther wire. This is easily achieved when L is large because 
the Josephson coupling strength increases proportional 
to the length of the coupling region, while the charging 
energies fall off as 

In Sec. lXIl we address the robustness of the topological 
degeneracy to local fermionic perturbations in the middle 
of the wire. Here, the situation becomes more compli- 
cated since back-scattering generates kinks in the super- 
conducting phases of various channels which are identical 
to impurity generated quantum phase-slips. In the case 
where the number of wires in the bundle is even, 
we find that backscattering, even in the center of the 
wire bundle, is capable of destroying the exponential in- 
dependence of the energy splitting. When the number of 
wires is odd, this does not happen, but any backscatter- 
ing that is present due to violations of momentum conser- 
vation near a junction between two wires will still lead to 
destruction of the exponential dependence of the energy 
spitting. (Fig. [^Jb)) Therefore, topological degeneracy 
that is exponential in the length of the wire can only oc- 
cur if all wires are in contact with a nanowire with a very 
large number of channels. The different responses of to- 
tal odd and even number of channel systems to quantum 
phase slips can be understood in terms of the fractional 
Josephson effect shown by TS wires with odd number of 
channels^^^^^ 

In Sec. lXIIl we find that there is another important dif- 
ference in the spectra for tunneling into the end of a TS 
wire which distinguishes the proximity-induced supercon- 
ductor from the isolated wire. In the proximity-induced 
case, there will be a finite amplitude for tunneling into 
the zero energy state, even in the limit of L — >■ 00. For 
the isolated wire, the amplitude for tunneling into the 
lowest energy state at the end of a wire falls off as a 
power of the length L. The missing spectral weight is di- 
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verted to states in which varying numbers of low energy 
"phonons" ( These phonons of course are the collective 
bosonic modes of the electron system, not of the underly- 
ing lattice, with a gapless linear energy dispersion in mo- 
mentum, arising from the short-range electron-electron 
interaction.) are simultaneously excited. In the limit 
L — >■ oo, this becomes a continuum with a spectral weight 
diverging at zero energy. 

Finally, in Sec. IXIIIl we show that despite the pres- 
ence of gapless phonons that broaden the end tunneling 
peaks associated with MFs, it is still possible to manip- 
ulate the topological degenerate space in wire bundles. 
This is done using a tunneling protocol which changes 
the tunneling between MFs in such a way that the tun- 
neling pins the local superconducting phase around the 
MFs before introducing tunneling between the MFs. 



III. GENERAL CONSEQUENCES OF NUMBER 
CONSERVATION 

To understand the difference between the isolated wire 
and the wire with proximity-induced superconductivity, 
recall that the two low-lying states formed from the MF 
modes at the wire ends correspond to states of opposite 
fermion parity for the wire. For the proximity-induced 
situation, tunneling of Cooper pairs will cause large fluc- 
tuations in the number of electrons in a long wire, while 
preserving the fermion parity. When the number fluctua- 
tions are large compared to unity, the difference in energy 
between an even and odd parity state will become expo- 
nentially small in the length of the wire. By contrast, 
for an isolated wire, where the Hamiltonian commutes 
with the total electron number, any energy eigenstate of 
the wire will be characterized by a definite particle num- 
ber N. For small deviations of N from some reference 
number Nq, we may write the ground state energy as 

Eo{N) = Eo{No) + a{N - No) + b{N - iVo)^ + ... (1) 

The constant b may be identified with the "charging en- 
ergy" for the wire; it is inversely proportional to the com- 
pressibility and it falls off as l/L, for a long wire with 
short-range interactions. The constant a depends on the 
reference energy for an isolated fermion, or on the chemi- 
cal potential of a reference reservoir, so it has no intrinsic 
meaning. The proper definition of the energy gap for tun- 
neling into the ground state of the wire is one-half the 
difference in the energies necessary to add or subtract 
an electron from the wire, which is simply equal to b in 
the present case. We note that the lowest energy for a 
particle-hole excitation in a finite wire will generally be 
the energy to excite a phonon in the mode with wave- 
length L/2, which also varies as l/L for large L, and 
is of similar magnitude to b. Therefore to achieve true 
topological degeneracy it will be necessary to move away 
from the paradigm of an isolated wire as will be done in 
SecHX] 



IV. NANOWIRE HAMILTONIAN AND 
BOSONIZATION 

While an isolated wire cannot show true topological 
degeneracy, it is instructive to study isolated nanowires 
before analyzing systems of wires. An isolated spin-orbit 
coupled multi-channel nanowire with Zeeman splitting 
has a band-structure shown in Fig. [TJ For a system with 
weak interactions, one can linearize the fermionic excita- 
tions in the vicinity of the Fermi energy fi 

^t(^) ^ e^'=---V'lLW + e-"=--^VU(^) (2) 

where fc^.a and a-re the fermi wave-vector and 

fermionic operator associated with band a — 1, . . . , Nch 
while 4'a ni^) ^-nd V'^ represent fermionic excitations 
around k « ±fci? a respectively22i2^. The linearized non- 
interacting part of the Hamiltonian is written as: 

Ho~^ dxil)\^,;^{x){iVF.adx)^a,L{x) 
a 

- fpi.R{x){iVF,adx)'lpa,R{x). (3) 

The attractive interaction between the fermions can be 
written as a product of 4-fermion operators 
and conserves total momentum so that only 2 classes of 
terms are allowed. The first of these classes 

Hi = J dx ^ Vry,a,bPaA^)Pb,r'ix), (4) 
a,fc,r,r' 

where pa,r{x) = tp\ j.{x)ipa,r{x), involves the densities of 
left and right movers in each band. The other class is a 
back-scattering term 

/ dxj2ylT'^i^lLi^)^lRi^)M^)M^). 

(5) 

which may be thought of as transferring Cooper pairs 
(V-'I ^(a;))from band b to a. The latter term will 

be shown to lead to a quasiparticle gap of the supercon- 
ducting kind when Hi contains attractive interactions. In 
this paper, we analyze the quasi-one dimensional Hamil- 
tonian H = H0 + H1+H2 using the Abelian bosonization 
identities^ii^^i^Si^ to represent the fermions in terms of 
bosonic fields as: 

i^aRix) = ^e-i^^"e'^^"(-) (6) 
V2A 

together with tpaLix) = -ipaR{.-x) where ipaR^x) = 
ipaRix + 2L) is assumed to be periodic and defined from 
[—L,L], L is the length of the wire and A is a short 
length-scale cut-off scale. The operator ANa is the num- 
ber fiuctuation operator that is associated with the band 
a defined as 

A^, =N^-^, (7) 
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where Na is the number of electrons in channel a, and 
the Klein factors F„ are constructed for this case as 



AN„ 



(8) 



to ensure canonical anti-commutation relations of the 
fermion operators^. The bosonic field -da (x) has a some- 
what involved commutation relation^^, which can be sim- 
plified by rearranging "daix) into a pair of fields 9a{x) and 
4>a{x) defined on the interval x S [0, L] by 



eaix) 



(t)a{x) 



^a{x)-'da{-x) 



(9) 

which are similar to the corresponding fields defined in 
bosonization of systems with periodic boundary condi- 
tions. The fields Oa{x) and (j)a{x) are periodic and anti- 
periodic respectively and have the commutation relations 

[eaix),ebiy)] - [Mx),My)] = o and [eaix),Mx')] = 

«7r[sgn(a; — x') — 1 -I- ^] so that na{x) = dx4'a{x), which 
is the density fiuctuation operator {pa,R{x) + pahix)) can 
be considered to be canonically conjugate to the variable 
9a{x) SO that 



[Tra{x),9b{y)] = iTr[S{x - y) - -r]Sab- 



(10) 



The 1 /L term arises from the absence of fc = modes of 
TTa{x) and 9b{x). The Cooper pair operator is given by 
i^iR{x)i^l^{,x) cx e^'''-^^^ so that 2ea(x) is related to the 
fiuctuating superconducting phase. The transformation 
of the Hamiltonian to a bosonic representation in terms 
of 9a{x) transforms the combination of term Hq -\- Hi in 
the Hamiltonian into an exactly solvable quadratic form 

Ho+Hi= JdxY, Uab^a^b + Vabdead9b (11) 

where Uab = VF,aSab + VR^R^a,b - Vn^L.a.b and Vab = 
VF.aSab + VR^R^ab + VR^L,a,b- The Coopcr pair scattering 
term H2 in the bosonic representation takes the form 



a^b 



dxV^ 



{pa 



a,b 



•cos{2{9a{x)-9bix))}, (12) 



which, for attractive interactions (V^'"^"*'^'' < 0), tends to 
align the phase variables 9a{x) between channels. Note 
that in realistic systems, the spin-orbit coupling is critical 
for generating the Cooper pair scattering term y^^j^^^'^— 
and the Zeeman splitting controls the number of channels 



V. TOPOLOGICAL SUPERFLUID ORDER 

In one-dimension, weak attractive interactions between 
fermions do not lead to a long-ranged superfiuid or- 
der parameter. However the Cooper pair tunneling 



term _ff2, becomes relevant under the renormalization 
group (RG) for attractive interactions Hi so that the 
phase differences 9a{x) — 9b{x) develop long-range or- 
der {9a{x) — 9b{x)) = or TT. To see this, we estimate 
the scaling dimension of the Cooper pair tunneling op- 
erator Pabix) = e'^^{^'^'y^)-^b(x)) ^ ^jj^^ jg contained in H2. 
The scaling dimension of Pabix) is estimated by calculat- 
ing the bulk correlator^ {Pab{x)P2b{v)) using periodic 
boundary conditions and the unperturbed Hamiltonian 
-ffo + Pti which is quadratic in the fields 9a{x). The cor- 
relators can thus be found using Gaussian integration 
to be {Pab{x)P:b{y)) = eMl dkj:^^^ CMk)9^k))Cp) 
where we have used (a = '2e^''^6a,a — '^e'^'^^Sb^a- Substi- 
tuting, and integrating over k we get {Pab{0)Pab(^)) ~ 
g-aC Cloga; ^ where we have used the matrix- 

vector notation — (.c ^ ^Sa^c — 26b^c and 



^T(^-l/2^^-l/2)-l/2^ 



< 1 



(13) 



provided that C^(y - U)( = 2C^Vr.lC < 0, i.e., the 
interaction Vr^l is attractive in all channels. Thus the 
Cooper-pair operator Pab{x) has a scaling dimension —a, 

with a < 1, so that V^^j^^^^ is relevant under RG and 
flows to large values at long length-scales. The phase 
differences 9a{x) — 9b{x) develop long-range order [i.e., 
{9a{x) - 9b{x)) = 7r(7ia - «&) for integers na,nb = 0, 1), 
so that 

9a{x)^na^ + 9{x)+59a{x) (14) 
where 59a{x) ^ tt over long length-scales. 

A. Ising pseudo-spin representation 

The addition of a single fermion -0^ (x) , which contains 
a factor e*'^"*-^-' (Eq. [5]), introduces a kink in the phase 
9a{x) so that 9a{x) — 9b{x) changes by tt. Such a kink 
costs a finite energy leading to a bulk quasiparticle gap 
in the superfiuid. However, there is no gap to tunnel- 
ing fermions at the ends since these do not introduce 
kinks. To describe fermion tunneling via such kinks it 
is necessary to accomodate a position dependent na{x) 
for a = 1, . . . ,Nch- Furthermore, the fermion tunneling 
operators contain fields 0a (x) and ANa which are canon- 
ically conjugate to na{x). 

Despite the bulk quasiparticle gap, the ground state 
of H is degenerate in the bosonic Hilbert space spanned 
by the operators 9a{x), since the Hamiltonian H is in- 
dependent of the value na = or 1 for each channel a. 
As discussed in the introduction, this is similar to the 
Z2 order parameter of the one-dimensional Ising model 
with a weak transverse field, where the system is char- 
acterized by an order parameter degeneracy (which is ria 
in this case) which takes 2 values. This correspondence 
is made more explicit by defining the local spin-fields for 
each channel 



S'f\x) 



(15) 
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The TS wire in this paper is in an ordered Ising phase 
{Si°'\x)) = ±1. Furthermore, the wire has a Z2 Ising de- 
generacy for each of the Nch channels. Furthermore, the 
low-energy effective Hamiltonian H (where spatial varia- 
tions S0a{x) have been integrated out), can be separated 

H[{S(^{x)},e{x)]^He[9{x)]+H,,[{Si'^\x)}] (16) 

where a — 1, ... , Nch and Hg[9(x)] arises from the com- 
mon phase variable 9{x) and keeps track of the global 
charging energy, while Hpg is the pseudo-spin (ps) Hamil- 
tonian keeping track of the degenerate sector described 
by the integers ria = or 1, which we will map to a 
pseudo-spin degree of freedom in Sec. I Villi The Hamil- 
tonian He[9{x)] describes massless phase fluctuations of 
9(x) and therefore corresponds to a low-energy effective 
action 



Zn 



dxdr 



{dr9{x,T)f 



+ u{d^9{x,T)f 



(17) 



where K is the Luttinger parameter associated with the 
phase variable 9{x) and u is the corresponding phonon 
velocity. 

The one-dimensional transverse field Ising model is re- 
lated to the mean-field BCS Hamiltonian of a topologi- 
cal spinless p-wave wire^ by a Jordan- Wigner transfor- 
mation. The order parameter degeneracy of the Ising 
model is mapped to the topological degeneracy of the 
p-wave wire by the non-local Jordan- Wigner transfor- 
mation. For our system, the bosonization transforma- 
tion can be thought of as a generalization of the non- 
local Jordan- Wigner transformation that allows us to de- 
scribe the topological degeneracy in terms of an order- 
parameter degeneracy. However, it is not a priori ob- 
vious that the order-parameter degeneracy associated 
with the ria integers in our system is robust to all lo- 
cal fermionic perturbations. We will establish that this 
is essentially true in Sec. IVIII and the following sections. 



VI. FIXED NUMBER HILBERT SPACE 

The physical Hilbert space of the nanowire system ex- 
tends only in a subspace of the bosonic Hilbert space 
with a fixed number of fermions Ntot- Moreover, the 
many-body wave- function is 27r periodic in terms of 9a (x) 
and not just 9{x). Thus, the many-body wave-function 

, when written in terms of the variables ria labeling 
the ground state sector and the average phase variable 
9 = J ^9{x), has the symmetry 



^[{na + l},9 + TT]='f[{na},i 



(18) 



which means that the wave function is unchanged if we 
change 9 by n and simultaneously change all the ria by 
1 for a = 1 , ... , Nch ■ Since the average phase 9 is canon- 
ically conjugate to the total fermion number Ntot, the 
wave-function has the form 



The symmetry condition Eq. [18] reduces to Ntot [{'^a + 
1}] = (-l)^*°*^'jv,„J{7ia}], where (-1)^"" is the fermion 
parity. Furthermore, since the operator (—1)^^" = 



has the effect of transforming ria 



1 (note 



[ANa,na] = i), the number conservation constraint has 
the simple form 



-1) 



(20) 



The Hilbert space constraint from Eq. [20] is simply that 
the physical states must be eigenstates of the operators 
ANa from the pseudo-spin space and Ntot with the 
same eigenvalue parity. 



VII. LOCAL END PERTURBATIONS 

As discussed previously, a degeneracy can be called 
topological only if it is insensitive to all local perturba- 
tions. So far, by bosonizing the Hamiltonian, we have 
identified an order parameter (i.e., the integers ria) in the 
bosonic space which may potentially lead to a topolog- 
ically protected degeneracy in the fermionic representa- 
tion. However, it is not a priori obvious that the bosonic 
order parameters ria are insensitive to any local pertur- 
bation in the fermionic representation. For example, im- 
purity scattering, which breaks momentum conservation, 
can introduce new terms in the Hamiltonian, such as 
ipl j.{x)'iph^r'ix). These terms are sensitive to the local 
value of the integer (ria — nb) associated with the ground 
state. The fermion operator -0^ introduces a kink in the 
phase variable 9a{x) that in general creates an orthog- 
onal excited state which is gapped by the kink energy. 
Therefore kinks near the ends of the wire that arise from 
end scattering are the dominant scattering terms that 
control the topological degeneracy in a wire with an odd 
number of channels. However, as we shall see in Sec. IXIl 
in a cluster of coupled odd-channel wires, if the number 
of coupled wires is even in some region of space, bulk 
scattering can lift the topological degeneracy in a subtle 
way, despite being disfavored by the superconducting en- 
ergy gap. We shall see that topological protection of the 
degeneracy requires having a wire with a large number 
of channels in the bundle of wires. 

In contrast to bulk scattering, inter-channel scattering 
at the ends can have a large effect in general. This is 
because the end fermion creation operators involve the 
4>a{x) at the ends which vanish and therefore do not cre- 
ate any kinks in the phase differences 9a{x)—9b{x). In the 
phase representation, the end scattering terms become 

H scattering (X - 0) iFaFb COS {9ai0) - 9b{0)) 

- iFaFbe'''^"-^-"'''> (21) 
at a; = and 

Hscatter^ngix ^ L) ^ iFaFbC'^^^'^^-^'''^ COS (0„(L) - 9b{L)) 

- iF„Fbe'''('^^"-'^^'')e"("''-"'') (22) 



*iVtoJ{"a},^] =e 



iNtate 



Nt. 



(19) 



at X 
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VIII. TOPOLOGICAL DEGENERACY 

The degeneracy associated with na = 0, 1 is in general 
lowered by Hscatteringix ~ 0,L). Only degeneracy that 
survives all weak local perturbations qualifies as being 
topologically protected. For simplicity of discussion, we 
will ignore the bulk scattering disorder, which we have 
shown does not affect wires with an odd number of chan- 
nels. As defined in Sec. IV Al the low-energy effective 
Hamiltonian in the degeneracy space spanned by is 
conveniently defined in terms of the pseudo-spin matri- 
ces 5*1"^ = e™"'^ ""'^ — aiTrANi 

sumed na{x) — 
of X in Eq. [TS] This is justified for the ground-state 
of a Hamiltonian of the form in Eq. [1^] where there is 
no term creating kinks in the order parameter na{x). 
The Klein-factors can be written in terms of spin op- 
erators as Fa = Y[b<a^^^^- '^^^ scattering terms in 
the spin- representation are written as Hscatteringix 



iriaTT g^^^ _ giTTAjVa^ where we have as- 
Ua to be free of kinks and independent 



0) 

S: 



(a) 



qia) r.(b) 



n 



'' na<c<6 ^■ 



a<c<6 



S^^\ and H, 



scattering 



[x ^ L) 



It follows from Eq. [20] that the to- 



tal fermion-parity operator can be written as 



(23) 



which commutes with all terms in the Hamiltonian. Thus 
the eigenstates of the Hamiltonian can be split into 2 sets 
of even {Q = 1) and odd [Q = —1) number of fermions. 
Topological degeneracy of even and odd fermions, simi- 
lar to that associated with MFs, occurs if there exists an 
operator A that commutes with all terms in the Hamil- 
tonian H and anti-commutes with the fermion parity op- 
erator Q. In this case 

|Q = -1) =A|Q = 1) (24) 

E{Q = -1) \Q = -l)=H\Q = -1) =HK\Q = 1) 

= KH\Q = l)= AE{Q = 1) IQ = 1) = E{Q = 1) |Q = -1) 

(25) 



proving that 



E{Q = 1) = E{Q = -1). 



(26) 



The operator A cannot exist in general in the even num- 
ber of bands case. To see this observe that if A commutes 
with both the end boundary terms Hscatteringix ^ 0, L), 
it must also commute with the product of these operators 
Slf^si'l For even number of bands, the total fermion- 
parity operator Q is a product of such operators. There- 
fore A commutes with Q as well and cannot anti-commute 
with Q. 

On the other hand, for odd number of bands, Q is a 
product of an odd number of 5*" and therefore A does 
not in general commute with Q. In fact, the operator A 
can be explicitly constructed as: 



A 



2b+l) 



(27) 



(a) 



a 





1 ■ ' ' 




— Interaction induced MF 

- ■ Proximity-induced MF 






, \ — 



(b) 



V (arb units) 



MF, V y MF4 



FIG. 2: (Color online) (a) Tunneling conductance G{V) into a 
Majorana site at the end of a long wire, shows power-law peak 
(solid line) characteristic of tunneling into Luttinger liquid 
with attractive interactions. This differs qualitatively from 
the (5-function peak (dashed line) expected from a proximity- 
coupled nanowire.— For the plot, we have chosen a Luttinger 
parameter K = 2 and added a lead-induced lorentzian broad- 
ening {{lu'^ -l-F^)"^) with width F = 0.02. (b) Schematic 
figure to obtain a topological degeneracy using topological su- 
perfluid wires containing phase fluctuations and phase-sfips. 
Three topological wires (shown in black) each of which sup- 
port MFs (discs at ends) are coupled by tunneling (shown as 
shaded regions) to obtain a topological degeneracy which is 
robust to bulk scattering and respects number conservation. 



which anticommutes with Q since a takes an odd num- 
ber of values. Also A commutes with both boundary 
scattering terms in H. Thus H has a topologically pro- 
tected degeneracy between even and odd number fermion 
modes only in wires with an odd number of channels. 
From here onwards we refer to such wires as topologi- 
cal wires. This proof is very similar to the proof of the 
Lieb-Shultz-Mattis theorem by Altman and Auerbach^. 
Moreover, this result is robust to bulk scattering (Eq. 1^5)) 
since this only effects even number of channels. How- 
ever, in the physical Hilbert space determined by the 
constraint Eq. 1201 which in the spin-representation re- 
quires that Q = (—1)^'°' where Ntot is the total number 
of fermions, the degeneracy sub-spaces Q — ±1 must 
have different energies because of Hq [Ntot] ■ This lifts the 
topological degeneracy even for an isolated wire with an 
odd number of channels. 



IX. TOPOLOGICAL DEGENERACY IN 
MULTIPLE WIRES 

The topological degeneracy in odd-channel wires is 
broken by the global phase fiuctuations which are needed 
to enforce number conservation. This makes an isolated 
wire unusable as a topological qubit. The number con- 
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servation constraint (Eq. [19)) can be relaxed by using a 
system of three odd channel topological wires which are 
connected in the middle by tunneling (Fig. [2lb)). The 
third wire is simply to keep the number of channels odd in 
the contact region so as to avoid bulk scattering (Eq. [28| . 
For simplicity of discussion we will focus on the other 
two wires with end MFs MFi_2.3,4- Considering specif- 
ically the case where the two wires have different fermi 
wave-vectors, so that single-particle tunneling between 
the wires is avoided, the dominant tunneling process be- 
tween the wires is Cooper pair tunneling. The overlap- 
ping segment of wire in the middle [-^ — Li, + Li] can 
be taken to have a common global phase d{x) so that 
one can define pseudo-spin variables ria for each wire 
j — 1,2. The geometry in Fig. [2][b) is such that there 
is no direct tunneling between MFs in the two different 
wires. Therefore even though the relative phases of the 
wires fluctuate away from the overlap region, the tunnel- 
ing Hamiltonian in each wire is unaffected and can be 

written as H = He[Ntot] + i?ps["i^^] + Hps[n'^^] where 
the constraint on Nt^t is now related to the total fermion 
parity Q = QiQ2- Moreover, the fermion parity opera- 
tors (5i,2 of individual wires are decoupled and the states 
with the same total fermion parity {Qi,Q2) = ("lil) 
and (Qi,Q2) = (I7— 1) are exactly degenerate up to ex- 
ponential splittings from instanton effects. 



X. EXPONENTIAL SPLITTING 

The pseudo-spin Hamiltonian Hps which leads to topo- 
logical degeneracy ignores virtual processes that over- 
come the tunnel barrier to relative phase fluctuations im- 
posed by the term V^l°'^^'^ J cos2{6'a(2:) — 9b{x)}. Such 
processes lead to an instanton tunneling"^^ between the 
states with different values of the local order-parameter 
na ~ njj — ±1 in the Ising representation, and can be 
treated as a tunneling of the fc = mode of 9a{x) — 6b{x) 
from one energy minimum to another'^^. The tunneling- 
induced splitting of the topological degeneracy can then 
be calculated by considering the Hamiltonian of the k = 
mode of the bosonized Hamiltonian H = Hq + Hi + H2 

and is found to scale as^S (5 - e'^'-^V^^^^^^. 

Similar splitting of the degeneracy occurs from phase 
slips between different wires which would occur in the 
limit of two proximity-induced TS nanowires, with a 
weak superconducting link between the two associated 
bulk superconductors. Again, we have energy splittings 
induced by coherent phase slips between the supercon- 
ductors, and the energy splitting should decrease expo- 
nentially with the square-root of the ratio between the 
Josephson coupling and the charging energy of the two 
superconductors. For bulk superconductors, the latter 
should be determined by the Coulomb charging energy, 
rather than the Fermi energy divided by the wire length 
or volume. 



XI. LOCAL POTENTIALS AND PHASE-SLIPS 

In Sec. IVH[ we argued that the effect of local poten- 
tials was dominated by scattering at the ends of the wire. 
Therefore, so far, we have established a topological de- 
generacy that is robust to local perturbations at the ends 
of the system. In this section, we examine more closely 
the effect of local potentials in the bulk of the wire. Such 
local potential terms of the form gilj'l{x)'4'b{x) introduce 
kinks in the order parameters Ua (x) , Uf, (x) , with coupling 
constant g, that in general leads to an excited state which 
is gapped by the kink energy. The kink generated by 
such a fermion tunneling term in the general case, is as- 
sociated with an excitation energy that takes us outside 
the ground-state manifold and therefore can only par- 
ticipate in virtual excitations that generate low-energy 
effective Hamiltonian terms which are products of tun- 
neling terms. The only such low-energy effective term 
that can be generated by bulk scattering is 

that generates a kink in all the channels a so that no rel- 
ative phase difference is created and the resulting wave- 
function can be described as a kink in the field 9{x). The 
term e*'*^^) in Eq.[28] is a vortex creation operator in the 
bosonic phase field 9{x). The corresponding term in the 
action, after the phase degree of freedom has been inte- 
grated out,^^ can be written as 

/„l/2 A K \ ^^1^ 

-(^) (-1)^"""^^^ (29) 

where the Luttinger parameter K defined in Eq. I17[ and 
K = K/Nch > 1 is a number of order unity. Such a term, 
being a product of Nch/'2, tunneling terms is exponen- 
tially small in the number of channels. Moreover, since 
each tunneling term has a pair of fermion operators, such 
a term can only be generated in a nanowire with an even 
number of channels. The bulk scattering term couples 
the pseudo-spin degree of freedom na{x) to the gapless 
phase modes 4>{x) making the even channel case useless 
for TQC. Such bulk scattering is avoided by introducing 
a wire with a very large number of channels, as shown in 
Fig. mb), so that the product term Eq. [51] is exponen- 
tially small. Thus, in order to avoid the bulk scattering, 
it is necessary to couple the nanowire to a superconductor 
with a large number of channels. This effectively requires 
the introduction of a bulk-like superconductor as in the 
case of the proximity effect. 

The bulk-scattering term in Eq.[2H]can be thought of as 
an impurity generated quantum phase-slip. The effects of 
a phase slip center can be understood, alternatively, by 
starting from the limit of strong backscattering, where 
the center effectively divides the wire bundle into two 
almost-separate halves. If the total number of channels 
in the wire bundle is odd, then dividing the bundle in two 
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pieces leads to a new Majorana state on each side of the 
spht, with a smah coupUng between the two Majoranas 
due to the residual tunneling. This coupling leads to a 
finite energy splitting between states where the resulting 
fermion state of the pair is occupied or not. If one intro- 
duces a change of 2tt in the phase of the superconductors 
on one side of the split, this will change the sign of the 
energy splitting, thus interchanging the upper and lower 
energy statesj ^'^^i^^'^^i^" If the phase slip occurs rapidly, 
the fermion number of the pair cannot change, so the 
energies before and after the phase slip will differ by the 
energy splitting between the two states. Thus, an iso- 
lated phase slip can only occur as a virtual process. (A 
double phase slip is allowed, however. Also, at finite tem- 
peratures, single 27r phase slips can occur, as a fermion 
at the center can be further excited into a bulk quasipar- 
ticle which is then eventually transferred to a low-energy 
state at the wire end.) If the Josephson coupling across 
the phase slip center is increased, the energy splitting 
between the occupied and unoccupied localized fermion 
state might become equal to the bulk energy gap, which 
only increases the action barrier to forming a single phase 
slip. 

For a wire with an odd number of channels, quantum 
tunneling of a single 27r phase slip is not allowed, and 
phase slips occur only in pairs, which have no net effect 
on the Majorana states. However, when an even num- 
ber of odd-channel wires is bound together to form a 
bundle with an even-number of channels, backscattering 
near the center of the bundle can lead to a 27r quantum 
phase slip. (Fig. EJb)) This has the effect of multiply- 
ing the Majorana operators to the left of the phase slip 
by a factor of -1 relative to operators on the right. The 
cumulative effect of such phase slips is to cause oscil- 
lations in any operator formed from a product of odd 
numbers of Majorana operators on opposite sides of the 
phase-slip center, and which signifies an energy splitting 
between various states in the Majorana fermion Hilbert 
space. [As an alternate description, we may say that the 
coherent generation of phase slips leads to a disruption 
of the continuous exchange of Cooper pairs across the 
phase-slip center, which leads in turn to a sensitivity of 
energy to the total electron parity of on one side of the 
center.] Similarly, backscattering near a junction of sev- 
eral odd-channeled wires can cause a 27r phase slip in one 
pair of wires, which would introduce a term in the energy 
sensitive to the total electron parity of the pair. Consis- 
tent with the bosonization results, we note that when 
the number of channels in each wire is large, the rate 
of generation of quantum phase slips will decrease expo- 
nentially with the number of channels, for a fixed small 
amplitude of single-particle backscattering, so that the 
energy splitting that is non-exponential in L will have 
an amplitude that exponentially small in Nch- In this 
way, we find that for a set of sufficiently thick wires with 
large but finite L, one can recover the exponential in L 
behavior found in the case of an adiabatic junction. 

We may ask further what happens if the tunnel cou- 



pling across a bulk phase slip center is taken all the way 
to zero. Suppose there is an odd number iV^, of odd- 
channel wires in the wire bundle, with each wire having 
a Majorana state at each of its free ends. Including the 
single Majorana state at the phase slip center on each half 
of the wire bundle, there will be an even number -f 1 
of Majorana states on each half of the system. Each wire 
will then have a number M = 2^^™"-'^)/^ of degenerate 
ground states, so the total system will have a degener- 
acy = 2^™~^. When couphng is reintroduced, this 
number stays the same. The pair of Majoranas near the 
phase slip center are removed from consideration as they 
now form a fermion state with finite energy. The re- 
maining 2N^, Majorana locations now belong to a single 
wire bundle, with one overall conserved charge, so the 
number of degenerate states is again 2^™~^. However, if 
the tunnel coupling across the center is very weak, the 
extra finite fermion state may have an energy which is 
well below the gap. If the total number of electrons in 
the system is specified, say, to be an even number, the 
energy splitting between the two fermion states may be 
thought of as arising from a difference between a state of 
the wire with a small excess electron probability near the 
center (and a small depletion near the ends) compared 
with one with a small depletion near the center. 

In the case where is even, there are no extra Majo- 
rana states near the phase slip center when the two halves 
are decoupled. The total number of zero-energy states is 
2N„-2 Tffi^Qj^ ii^Q wires are completely decoupled, while it 
would be 2^™~^ if there were no backscattering and no 
phase slips allowed. For the case of weak tunneling there 
would be two sets of states with 2^™"^ degenerate states 
each, but a small energy splitting between the two sets. 
This energy splitting corresponds to a splitting between 
states where the fermion parity on, say, the right hand 
side is predominantly odd and where it is predominantly 
even. (The fermion parity on the other side is accord- 
ingly determined since we are considering a total system 
with a definite number of electrons.) 

XII. END MFS AND TUNNELING SPECTRUM 

Even though the exact topological degeneracy may be 
broken by number conservation for single wires, the de- 
generacy for odd-number-of-channel wires is only split to 
order and it should continue to manifest itself as 
a near-zero energy peak in end-tunneling conductance, 
as has been predicted in previous work. To understand 
the spatial dependence of the tunneling, it is convenient 
to generalize our definition of the pseudo-spin used in 

Sec. IVIIII so that si°'\x) and Sy°'\x) are position de- 
pendent field operators. The order parameter operator 

si'^^ (x) was already defined as a position dependent field 
in Eq. 1151 The generalized x spin operators are defined as 

Si^^x) - e'^'l'a(x)+AN^x^/L] ^gjjjg ^jjg ^^^^ ^j^j^^ ^J^) ^ 

is canonically conjugate to Oa{x), and therefore 
generates a kink in na{x). The corresponding y pseudo- 
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spin operator is defined as 5'^°' (a;) = iSz {x)Sx (x). Us- 
ing the bosonization identity Eq. [S] in conjunction with 
the Eq. [T31 the fermion tunnehng operator at the ends 
can be decomposed as ipl{x ~ 0) ~ e^^^^'>^a,R{x) into 
a manifestly Hermitean (and therefore Majorana part) 
7a _r(-^ ~ i) ^ FaSy°'\x) and a phase fluctuation opera- 
tor e'^'^-*. Similarly, the left fermion tunneling operator 
ipai^ ~ i) ^ e'^'^^^7a.L(a;) can be shown to contain a 
Majorana part 'jI i^{x ^ 0) ^ FaSz'^\x). The opera- 
tors Si^^x) and Si^^x) are localized at the ends of the 

wire since the bulk of the si°'\x) ordered Ising phases 
are gapped. Moreover, the Hamiltonian Hps[na] in the 
pseudo-spin representation and the corresponding eigen- 
states \Q = ±1) are both real. Therefore in this eigen- 
basis, the left and right Majorana operators are given by 
7i = (Tx and = Uy respectively, where (Jx,y &rG the 
Pauli matrices in the pseudo-spin basis where \Q = ±1) 
are eigenstates of CTz. While tunneling into these ideal- 
ized MF modes would have generated a sharp zero-energy 
peak similar to the proximity-induced casSii^, the end 
tunneling conductance into the physical fermions is domi- 
nated by the phase mode d{x) and is essentially identical 
to tunneling into the end of a single spin-less fermion 
channel with attractive interactions^i^^: 

G{V) oc (30) 

where K > 1, so that there is a power-law divergence at 
y — >■ (Fig.[2ja)), instead of the 5{V) from the proximity 
effect theory without phase fluctuations. 

More interestingly, for a finite wire of length L, the 
phonon spectrum is gapped with energy levels given by 
LOn = {n + l)hu/ L, where u is the phonon velocity de- 
fined in Eq. 1171 Such a finite wire leads to a discretized 
spectral function (instead of the continuum one shown in 
Fig. HJa)) consisting of a sequence of narrow peaks re- 
fiecting the phonon frequencies. One can also consider 
the possibility of tunneling a single fermion into the Ma- 
jorana end state from a neighboring confining potential 
such as a quantum dot^. Coupling the Majorana end 
state to such a localized fermion state not only switches 
the fermion parity Q but also excites a wave-packet of 
phonons that propagates from the end. Such a phonon 
packet can only be avoided by performing the tunneling 
operation on a time-scale that is adiabatic compared to 
the inverse lowest phonon frequency loq^ ^ L/hc which 
scales with the length of the system. 



XIII. TOPOLOGICALLY PROTECTED 
OPERATIONS 

In earlier sections we have shown that a pair of topo- 
logical wires have a topological degeneracy which is only 
split by an amount that decreases rapidly with the length 
of the wire L and therefore may be used as a topologi- 
cally protected qubit. However, the question of whether 



this topologically protected qubit can be used to perform 
topological quantum computation is not obvious. In par- 
ticular, the presence of phase-fluctuations and phase-slips 
can potentially lead to decoherence of otherwise topolog- 
ically protected manipulations of the qubit. Here, we 
will restrict our discussion to qubit read-out and MF 
braid operations, which are the 2 classes of topologi- 
cally protected manipulations that have been proposed 
for topological systems with Majorana fermions. While 
not sufficient to perform topological quantum computa- 
tion by themselves, they are enough to verify directly 
the non-Abelian statistics of MFs. Recently several 
schemes have been proposed for performing the opera- 
tions of read-out of topological qubits2ii2ii2^ and braiding 
of MFs.-l"^'' The proposed read-out schemes that use ei- 
ther interferometry-based fermion parity measurement 
capacitive charge read-out , or tunneling.— On the 
other hand, the braiding schemes require creation and 
annihilation of tri-junctions in the wires. All the read- 
out operations and the creation and annihilation of tri- 
junctions can in principle be implemented for the topo- 
logical superfluid wire systems discussed in this paper 
by a straightforward generalization of the correspond- 
ing proposals. However, some care must be taken when 
implementing these proposals to avoid effects due to 
phase-fluctuations and phase-slips during the read-out 
and braid operations, which can lead to unwanted de- 
coherence. The additional requirements, together with 
effects of finite phonon propagation velocities, can lead 
to limitations on the rate at which manipulations can be 
performed. A full discussion of these issues would re- 
quire a long digression, which we postpone to a future 
publication. 



XIV. CONCLUSION 

In this paper we have shown how a topological degener- 
acy emerges in a number conserving Hamiltonian where 
superfluidity is induced in a multi-channel nanowire with 
an odd number of channels by attractive interactions. 
Our treatment uses the bosonization formulation and 
therefore goes beyond the BCS mean-fleld -I- fluctuation 
treatment in a systematic way so that number conserva- 
tion is clearly preserved. In the bosonization representa- 
tion, which is related to the fermions in a non-local way, 
the topological order of the fermionic Hamiltonian finds 
a natural interpretation as a Z2 pseudo-spin order of the 
relative phases between various channels. However, the 
apparent degeneracy of 2^'='', where Nch is the number 
of channels, is explicitly shown to be broken down by 
end scattering to 2^<:''(™°''^). Therefore we show in the 
general interacting case that the topological degeneracy 
survives end scattering only for wires with an odd num- 
ber of channels. Local perturbations in the middle of the 
wire introduce quantum phase-slips that are shown to 
lift the topological degeneracy in general by an energy- 
splitting that scales as g^<:h/2j^-K.N^h ^ where g measures 
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the strength of the impurity potential and k is a constant 
of order unity. Thus the topological degeneracy which is 
robust to all local perturbation occurs only for a bundle 
of wires where one of the wires has a large number of 
channels. 

The number conservation constraint appears in a sim- 
ple way as a constraint between the pseudo-spin and 
phase sectors of the Hamiltonian. This constraint in 
general lifts the topological degeneracy of the topolog- 
ical nanowires unless multiple nanowires are considered. 
Despite the 1/L splitting of the topological degeneracy, 
an isolated odd channeled nanowire shows a peak in the 
tunneling spectrum. Unlike the proximity-induced super- 
conducting case, the peak is strongly broadened by phase 
fluctuations. Despite the presence of phase fluctuations, 
we find nearly perfect degeneracy for the Majorana states 
in a multiwire bundle, provided that phase slips are suf- 
ficiently suppressed in wire sections containing an even 
number of channels. 

As expected, our results approach the mean-field BCS 
limit in the case where the topological superfluid wire is 
in contact with an infinite bulk superconductor, which 
creates an ensemble of number states in the wire. In this 
case, a system of wires which has n MFs has a degen- 
eracy of 2"/^. The bosonization approach in this paper 
allows us to study how the limit of an infinite system 
is approached and explicitly consider the role of global 
and local number conservation, which are not manifestly 
obeyed by the mean-field formalism. Imposing global 
number conservation on the mean-field results splits the 
degenerate ground-state manifold into 2 groups of 2"/^"^ 
degenerate states with definite number of fermions in 
each state. The splitting between the 2 groups of states 
vanishes as an inverse power of the system size in the limit 
of an large system. Therefore, the splitting in the topo- 
logically degenerate manifold from global number con- 
servation vanishes in the limit of an infinite system in 



agreement with the BCS mean-field limit. The power 
law dependence on system size is too slow to be useful 
for quantum computation. This is not a problem, how- 
ever, as manipulations such as braiding of Majorana sites 
can be performed within a 2^/^^^ dimensional subspace 
which has definite particle number. 

The power-law dependence of energy splittings due to 
global number conservation can be a serious obstacle to 
topological computation if at intermediate times the sys- 
tem is divided into isolated wire segments with number 
conservation for the individual segments. The resulting 
multiple energy splittings can then lead to rapid dephas- 
ing of states in the computational Hilbert space. 

Local number conservation, which is present in our 
superfluid Hamilto- nian, manifests itself in more sub- 
tle effects such as local phase fluctuations and quantum 
phase slips. While local phase fluctuations do not lead to 
any additional split- ting of the degeneracy, as discussed 
in Sec. IIXI phase-slips split the topological degeneracy. 
When the number of channels Nch becomes large, how- 
ever, these splittings become rapidly small, and they are 
typically much smaller than the spittings due to global 
number conservation. For fixed system length, we find 
that the splittings decrease exponentially with Nch, while 
for fixed large Nchi the splittings decrease as a large in- 
verse power of the system length. Therefore, our result 
approaches the mean-field BCS result in the limit of an 
infinite system size and an infinite number of channels, 
both of which can be realized in the proximity-induced 
superconductivity case from a large bulk superconductor 
as in the previous proposals J^^i^ 
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